Charge pumping in carbon nanotube quantum dots 
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We investigate charge pumping in carbon nanotube quantum dots driven by the electric field of 
a surface acoustic wave. We find that at small driving amplitudes, the pumped current reverses 
polarity as the conductance is tuned through a Coulomb blockade peak using a gate electrode. We 
study the behavior as a function of wave amplitude, frequency and direction and develop a model 
in which our results can be understood as resulting from adiabatic charge redistribution between 
the leads and quantum dots on the nanotube. 

PACS numbers: 73.63.Kv, 73.23.Hk, 73.63.Fg, 85.35.Kt, 72.50.+b 



The dynamics of charge transport at the level of single 
electrons or electron spins in nanoscale systems such as 
quantum dots is important from both applied and funda- 
mental viewpoints. A sensitive probe of the response of 
a quantum dot to changing external parameters (e.g. a 
gate voltage or magnetic field) is the dc charge pumping 
current generated in the absence of an applied bias [l| . A 
general result for open, non-interacting electron systems 
is the Brouwer formula which relates the pumped cur- 
rent to the scattering matrix of the system 0] ■ In inter- 
acting systems, correlations add complexity and modify 
the predictions. Recent theoretical work has considered 
pumping in quantum dots with weak interactions 0, 0] , 
in the Kondo regime [f| , in the Coulomb blockade regime 
[g[ , and for superconducting leads @, [1] • For the case of 
interacting electrons in one dimension (the Luttinger liq- 
uid), current quantization [9l. 1 1 0t| and pure spin currents 
[ll| have been predicted. 

Experimentally, the rich variety of predicted phenom- 
ena has remained largely unexplored [12l . [l3j . A promis- 
ing system to test the various theories of charge pumping 
is carbon nanotubcs, since all of the transport regimes 
mentioned above have already been realized in conven- 
tional nanotube devices 0, UH, EH- In this Letter we 
pursue the potential of nanotubes, demonstrating charge 
pumping in nanotube quantum dots up to frequencies 
of 2.6 GHz and present a theoretical description of the 
observed effects. 

While charge pumping could in principle be realized by 
modulating side or to p g ates that are capacitively cou- 
pled to the nanotubes [171 ] , we adopt a different approach 
in which a surface acoustic wave (SAW) pumps charge 
through a carbon nanotube [3 01 IMS!!. An individual 
nanotube is contacted on a piezoelectric quartz substrate 
by palladium source and drain electrode that are sepa- 
rated by 5 fiva (Fig. 1(c) lower right inset). A side gate 
electrode is used to vary the electrostatic potential in the 
nanotube. Several mm beyond each contact are two SAW 
transducers (Fig. 1(c) left inset) having resonant frequen- 
cies /saw of about 2.6 GHz and 544 MHz. The electric 



field accompanying the SAW pumps charge through the 
nanotube. While the use of SAWs restricts experiments 
to the transducers' resonant frequencies, it avoids direct 
capacitive coupling between the high-frequency and de- 
vice electrodes and enables a clear distinction between 
the signal due to the SAW (possible only at the trans- 
ducers' resonant frequencies) and rectified currents from 
radiated fields (possible at all frequencies). Details of 
sample fabrication, transducer operation and propagat- 
ing SAW fields were described in Ref . [l9j] . 

Fig. 1(a) shows the dc transport properties of the nan- 
otube at temperature T = 5 K. The approximate peri- 
odicity and large charging energy Uc ~ 10—15 meV 
observed indicate that the nanotube is divided into two 
(or more) sections such that the conductance is domi- 
nated by Coulomb blockade of a single electron puddle 
that is much smaller than the 5 /zm source-drain sep- 
aration. This also follows from the suppression of the 
linear-response conductance and additional features (ar- 
rows) observed in the differential conductance (see inset 
to Fig. 1(a)) [22j]. Fig. 1(c) shows the induced current, 
/saw, in the presence of a SAW field at 2607 MHz in 
the absence of a source-drain bias voltage. For a SAW 
velocity vsaw ~ 3200 m/s on quartz this corresponds to 
a SAW wavelength Asaw = ^saw//saw ~ 1.2 /jm fl9l |. 
When a low power Psaw is applied to the SAW trans- 
ducers, a dc current is induced whose direction alternates 
as a sensitive function of gate voltage V g . The peak-and- 
dip features in the current are clearly correlated with the 
Coulomb blockade peaks, and the current changes polar- 
ity on reversal of the SAW direction. 

These features were studied in more detail as a function 
of SAW amplitude and for both available SAW frequen- 
cies. Figure 2(a) shows /saw as a function of Psaw- 
The peak-and-dip features increase in magnitude and 
move outwards as Psaw is increased. When we plot the 

derivative of the current with respect to the SAW ampli- 
1/2 

tude Vsaw oc Psaw as m Fig- 2(c), it becomes apparent 
that the dominant features move essentially linearly with 
Vsaw- At the highest applied SAW powers (up to 20 dBm 



2 



-10 10 



1.05 -1.00 -0.95 




(a) 



-0.2 0.0 0.2 



FIG. 1: (a) Color scale representation of the current as a 
function of Vsd and V g . Inset: differential conductance around 
a Coulomb peak (dark = more conductive) (b) Current trace 
at V s d = —2.5 mV showing Coulomb blockade oscillations, 
(c) SAW induced current as a function of V g in the absence 
of an applied bias. In red is the current produced by the left- 
moving SAW at 2602 MHz with P SA w = -15 dBm. In blue is 
the current produced by the right-moving SAW at 2607 MHz 
with Psaw = — 10 dBm. Right inset: simplified schematic of 
the device. Left inset: photograph of one of the transducers, 
which is able to generate SAWs at frequencies of both 2607 
MHz and 544 MHz. Top inset: frequency dependence of the 
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or 100 mW) the features merge into a background of cur- 
rent in the direction of the SAW across the whole gate 
voltage range [HI]. 

Assuming that the dominant action of the SAW is to 
modulate the electrostatic potential at the nanotube, we 
can estimate the actual value of Vsaw from the linear 
evolution of the current features with amplitude and the 
gate efficiency deduced from the dc Coulomb blockade 
pattern. We find a SAW amplitude Vsaw — 150 mV at 
the maximum power Psaw = 20 dBm applied to the left 
transducer. This agrees with estimates from transducer 
impedance calculations and direct transducer transmit- 
tance measurements The inset to Fig. 1(c) shows 

Psaw as a function of frequency. A peak is observed at 
the transducer resonant frequency only, demonstrating 
that the observed current features do not result from di- 
rectly radiated fields or from photon-assisted tunneling 
(see also Refs. 0[jj|). 

All the observed behavior can be explained by a model 
in which disorder in the nanotube causes it to contain two 
or more localized electron puddles in series. The SAW 
modulates the energies of the states in these puddles and 
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FIG. 2: (a) Color scale representation of the SAW-induced 
current as function of V a d and V g , showing sign reversal as V g 
crosses a Coulomb blockade peak. The line trace is taken at 
Psaw = —8 dBm. (b) Calculated current using the model ex- 
plained in the main text, (c) Derivative of the SAW induced 
current with respect to the SAW amplitude as a function of 
V g and (normalized) Vsaw- The maximum SAW amplitude 
is estimated to be ~ 100 — 150 mV. Inset: dl/dVsAw as cal- 
culated in the model, showing the same linear peak splitting. 



adiabatic charge redistribution between the puddles and 
the leads gives rise to a dc current. We first discuss a min- 
imal model in which there are two puddles, each acting as 
a quantum dot with a single spinless level (see Fig. 3(a)), 
before introducing asymmetries in the tunnel couplings 
and multiple levels in the dots to better reflect the likely 
experimental situation. Such a simple two-level system 
provides a useful approximation to the Coulomb blockade 
physics of tunnel-coupled double dots; see e.g Ref. [23|. 
The only parameters that enter the model are the tunnel 
couplings r l and Tr of the left and right dots to left and 
right leads respectively, Vsaw, the phase difference <j> of 
the SAW between the two dots, and the tunnel coupling 
A between them. The instantaneous effective Hamilto- 
nian of this double-dot system can be written as 
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ei - iT L /2 A/2 
A/2 e 2 - iT R /2 



(1) 



For any periodic time dependence of the dot energies 
El and 62, the adiabatic current 0, H3, [H, [26[ from 
left to right can be calculated as a surface integral 0], 
I = ef J A P(ei,e 2 ) de\ds2, over the area A enclosed 
parametrically by the pumping trajectory in the (ei, £2)- 
plane. 

The response function R can be obtained from Eq. 1, 
see e.g. Ref. [26| , and in the zero- temperature limit takes 
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a straightforward form, R{e\,E2) = 

-(32/TT)A 2 T L T R ,(e 1 T R , + e 2 T L ) 

[2T L T R A 2 + 4e 2 r| + \e\T\ + (A 2 - 4 £l e 2 ) 2 + r 2 rfj 2 ' 

The response function is shown graphically in Fig. 3(b,c) 
for different parameter values. It exhibits one pro- 
nounced minimum (blue) and one maximum (red). For 
the pumping trajectory, we assume the SAW periodically 
modulates the levels according to 
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Here a,\ and a 2 are the coupling efficiencies of the gate to 
the two dots and <5 is a level offset parameter. The effect 
of the SAW on the tunnel barriers is thought to be negli- 
gible for the small SAW amplitudes considered here. For 
the case of symmetric coupling Tl = Tr (Fig. 3(b)), we 
see that the peaks of negative and positive pumping cur- 
rent, as seen in Fig. 2(b) correspond to elliptical trajecto- 
ries around the two triple points in the stability diagram 
of the double dot [H, [27| • In obvious notation we may 
represent these as the (0,0) ->■ (0,1) -> (1,0) -> (0,0) 
electron and (1, 1) — * (1,0) — * (0, 1) — ► (1, 1) hole cycles, 
respectively. 

One property of the pumping current that follows im- 
mediately is that it changes polarity when the direc- 
tion of the pumping contour reverses, i.e., when the 
SAW direction is reversed. Another is that for small 
SAW amplitudes, I oc V^ 2 AW sin </> @], because in this 
limit R(si,e 2 ) is approximately constant within the area 
A = 7re 2 Vg 2 AW sin0 enclosed by the trajectories. 

As V g is varied, the center of the pumping trajectory 
follows a straight line which for the simple case ct\ = a 2 
and 8 = is the solid diagonal in Fig. 3(b)-(d). In this 
case, the current traces I(V g ) exhibit symmetric sign- 
reversing features. If the dots are single-level and com- 
pletely symmetric, i.e., Fl = Fr and 6 = 0, then the cur- 
rent traces for different /saw are dependent only on the 
ratio A/r. This simple situation already matches the ex- 
perimental data quite well, as shown in Fig. 2(b), taking 
A/r = 0.5. Figure 3 also illustrates the effects of asym- 
metries. A nonzero level offset 8 shifts the diagonal line 
[see e.g. the dotted line in Fig. 3(b)], while asymmetric 
coupling to the gate, a\ ^ a 2 , changes the slope (dashed 
line). The result is an asymmetry between the positive 
and negative current peaks which is most pronounced in 
the weak pumping regime. For larger SAW powers, such 
that eVsAW — 6, the difference between positive and neg- 
ative peaks begins to average out. In the limit of strong 
pumping, the pumping trajectories encompass both pos- 
itive and negative parts of R(e\,e 2 ) whose contributions 
to the pumping tend to cancel. This explains the ex- 
perimental observation (Fig. 2(c)) that the opposite-sign 
peak pairs in the pumping current vs V g move apart lin- 
early in Vsaw : the peak current occurs where the integral 
of R is maximal as a function of eg, arL d this occurs at 
|e | ~ eVsAw 
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FIG. 3: (a) Schematic of a double dot in which the single- 
electron states are periodically modulated by the SAW with 
phase difference <j>. (b) Color scale representation of the func- 
tion R(ei,£,i) for symmetric coupling to the leads. The SAW- 
induced current is proportional to the integral of this function 
over the area traversed in [£1,63] space. For <f) = n/2 the tra- 
jectories are circles with diameters proportional to Vsaw; for 
<f> = tv/10 they are narrow ellipses. The ordered pairs (n,m) 
indicate the electron occupancy of the quantum dots, (c) 
Same as panel (b) for asymmetric coupling to the leads. The 
direction £0 corresponds to the case of equal coupling of the 
dots to the gate, ct\ = a^. (d) Calculation of R(e\,e-i) with 3 
levels on the left dot taken to have spacing equal to 2A and 
identical couplings to the level in the other dot. 



The qualitative behavior of the model is quite insensi- 
tive to details such as asymmetry in the tunnel couplings 
or a multiplicity of levels in the dots. In particular the 
experimental situation in which the conductance is domi- 
nated by a single quantum dot can be modelled by taking 
a single level on the right dot and multiple levels on the 
left, as shown in Fig. 3(d). The sign-reversing nature of 
the IsAW-Vg traces is still a robust feature although sig- 
natures of extra levels may appear for certain pumping 
trajectories [281 ]. 

Finally, we examine whether the experimental effects 
of changing Asaw are consistent with the model. Fig. 
4(a) shows /saw at /saw = 2607 MHz in the vicinity of a 
Coulomb peak. A slight asymmetry between positive and 
negative peaks is seen for small Psaw- In the model this 
is reproduced (see Fig. 4(b)) by, for instance, assuming 
an asymmetric coupling to the leads and a small offset 
6, corresponding to the situation depicted in Fig. 3(c). 
Other kinds of asymmetry such as ol\ ^ a 2 have a sim- 
ilar result. Fig. 4(c) shows the current at /saw = 544 
MHz, using the second pair of transducers with larger 
finger spacing. The current is now nearly always of neg- 
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FIG. 4: (a) Color scale representation of the SAW-induced 
current as a function of SAW power and V g for /saw = 2607 
MHz. (b) Calculated SAW current for a double QD model 
with an asymmetry between the tunnel couplings Vl and Tr, 
as discussed in the text. The parameters are Tl '■ T_r : A : 
S — 1.0 : 0.1 : 0.5 : 0.3, and the phase difference cj> = 7r/2. 
(c) SAW-induced current for the same resonance peak as in 
panel (a) for /saw = 544 MHz. (d) Calculated SAW current 
as in panel (b) but for a phase difference <f> — 7r/10 [2S| 



ative polarity. Fig. 4(d) shows the corresponding model 
prediction, which agrees very well. Here, the effect of de- 
creasing /saw is just to reduce the phase difference be- 
tween the dots (0 — > when Asaw becomes much larger 
than the device dimensions). This reduces the width and 
area of the elliptical pumping trajectory. The result can 
be understood qualitatively from Fig. 3(c) noting that 
the level offset 5 shifts the line followed by the ellipse 
center towards the negative (blue) regions of R{s\,e2)- 
For large SAW powers and circular pumping trajecto- 
ries, the enclosed area includes a domain where R(si, £2) 
is positive. However, after decreasing Psaw (smaller ra- 
dius) or /saw (narrower ellipse) the integral only includes 
regions of negative i2(ei,£2)i enhancing the asymmetry. 
Note that the enhancement of the asymmetry in /saw 
when decreasing SAW power or frequency is expected 
whenever the center of the pumping trajectory follows a 
line in [ex, £2] parameter space that is asymmetric with 
respect to the minima and maxima of R(si, £2), irrespec- 
tive of the precise form of this function. The results in 
Fig. 4 are therefore consistent with the model also when 
in the experiment i?(ei, £2) is significantly more complex 
than the model situations shown in Fig. 3. To completely 
describe the experiment (the model actually underesti- 
mates the amount of pumped current and does not ac- 
count for the fine structure observed at 544 MHz) would 
require a better knowledge of the microscopic details of 
the device than we have at present. 

We thank A. Aharony, O. Entin-Wohlman, and L. Lev- 
itov for discussions. MB has been supported by the UK 
QIP IRC (GR/S82176/01). VK has been supported by 
Latvian Council of Science, European Social Fund and 
German-Israeli Project Cooperation (DIP). 



[15] 



D.J. Thouless, Phys. Rev. B 27, 6083 (1983). [16 
P.W. Brouwer, Phys. Rev. B 58, R10135 (1998). [17 
I.L. Aleiner and A.V. Andreev, Phys. Rev. Lett. 81, 1286 
(1998). [18 
P.W. Brouwer, A. Lamacraft, and K. Flensberg, Phys. [19 
Rev. B 72, 075316 (2005). 
T. Aono, Phys. Rev. Lett. 93, 116601 (2004). [20 
J. Splettstoesser, M. Governale, J. Konig, and R. Fazio, 
Phys. Rev. Lett. 95, 246803 (2005). [21 
M. Blaauboer, Phys. Rev. B 65, 235318 (2002). [22 
J. Splettstoesser, M. Governale, J. Konig, F. Taddei, and 
R. Fazio, Phys. Rev. B 75, 235302 (2007). 
V.I. Talyanskii, D.S. Novikov, B.D. Simons, and L.S. 
Levitov, Phys. Rev. Lett. 87, 276802 (2001). [23 
D.S. Novikov, Phys. Rev. Lett. 95, 066401 (2005). [24 
R. Citro, N. Andrei, and Q. Niu, Phys. Rev. B 68, 165312 
(2003) [25 
H. Pothier, P. Lafarge, C. Urbina, D. Esteve, and M.H. 
Devoret, Europhys. Lett. 17, 249 (1992). [26 
M. Switkes et al, Science 283, 1905 (1999). 
J. Nygard, D.H. Cobden, and P.E. Lindelof, Nature 408, [27 
342 (2000). [28 
M. Bockrath et al, Nature 397, 598 (1999). 



M.R. Buitelaar et al, Phys. Rev. Lett. 91, 057005 (2003). 
Y. Wei, J. Wang, H. Guo, and C. Roland, Phys. Rev. B 
64, 115321 (2001). 

P.J. Leek et al, Phys. Rev. Lett. 95, 256802 (2005). 
M.R. Buitelaar et al, Semicond. Sci. Technol. 21, S69 
(2006). 

J. Ebbecke, C.J. Strobl, and A. Wixforth, Phys. Rev. B 
70, 233401 (2004). 

Y-S. Shin et al, Phys. Rev. B 74, 195415 (2006). 
The features in the dl/dV are tentitatively interpreted 
as signatures of the charging energy of a second (smaller) 
quantum dot in series, see e.g. J. Park and PL. McEuen, 
Appl. Phys. Lett. 79, 1361 (2001). 

W.G. van der Wiel et al, Rev. Mod. Phys. 75, 1 (2003). 
M. Biittiker, H. Thomas, and A. Pretre, Z. Phys. B 94, 
133 (1994). 

O. Entin-Wohlman, A. Aharony, and Y. Levinson, Phys. 
Rev. B 65, 195411 (2002). 

V. Kashcheyevs, A. Aharony, and O. Entin-Wohlman, 
Phys. Rev. B 69, 195301 (2004). 

W.J.M. Naber et al, Phys. Rev. Lett. 96, 136807 (2006). 
This also holds when the level spacing on the left dot is 
small compared to the temperature or Tl and the fea- 



5 



tures in R(ei,£2) merge into two horizontal bands of op- 
posite polarity. 

[29] In the model we assume a reduction of by a factor of 
five, corresponding to the fivefold increase of Asaw- The 



correspondence between <j> and Asaw is uniquely defined 
only when Asaw is larger than twice the distance between 
the centers of the two dots. 



